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ANISOPLANATISM  IN  ADAPTIVE  OPTICS 

Feng  Yuezhong,  Gong  Zhiben,  and  Song  Zhengfang 

Anhui  Institute  of  Optica  and  Fine  Mechanics  (Hefei), 

Chinese  Academy  of  Sciences 

Abstract:  The  paper  discusses  the  influence  of  the 
compensatory  effect  due  to  the  anisoplanatic  field  in  adaptive 
optics,  and  proposes  a  new  method,  a  phase  gradient  method  to 
reduce  the  error  in  anisoplanatism.  As  shown  by  analytical 
results,  this  method  can  significantly  increase  the  range  of 
angles  for  effective  calibration;  thus,  the  influence  on  the 
anisoplanatic  error  is  greatly  reduced.  The  authors  anticipate 
that  the  phase  gradient  method  can  be  expected  to  solve  beacon 
anisoplanatic  problem  in  adaptive  optics. 

Key  words:  adaptive  optics,  isoplanatic  angle,  phase 
gradient  method. 

I.  Introduction 

Uhen  calibrating  the  effect  of  turbulence  in  the  atmosphere 
by  using  adaptive  optics,  in  the  conventional  work  approach 
measure  the  wavefront  error  information  carried  by  the  light  beam 

1 


r 


/ 


coming  from  the  target,  and  apply  the  conjugate  calibration 
amount  to  the  emitted  light  beam;  thus,  the  phase  distortion  in 
the  channel  is  precisely  cancelled  out  when  the  emitted  light 
beam  traverses  the  same  light  path.  The  light  beam  providing  the 
wavefront  error  information  is  called  the  beacon  light  beam, 
which  comes  from  the  target  or  other  independent  light  sources. 
Related  to  the  beacon,  an  important  parameter  is  the  isoplanatic 
angle,  which  forms  a  certain  conical  region  in  the  atmosphere. 

In  the  direction  of  the  isoplanatic  angle  in  the  region, 
turbulence  induces  phase  fluctuations;  these  are  correlated. 

When  the  beacon  conforms  to  the  target  direction,  the  wavefront 
distortion  of  the  beacon  light  beam  will  fully  represent  the 
phase  error  along  the  transmission  channel.  When  the  beacon 
light  beam  and  the  target  light  beam  form  an  included  angle  not 
exceeding  the  isoplanatic  angle,  the  beacon  light  beam  still  can 
provide  precise  phase  information  for  the  transmission  channel. 
However,  when  the  included  angle  exceeds  the  isoplanatic  angle, 
the  turbulent  regions  traversed  by  the  beacon  light  beam  and  the 
emitted  light  beam  will  be  quite  different;  the  phase  information 
of  the  beacon  light  beam  cannot  completely  represent  the  phase 
distortion  along  the  emission  path.  Thus,  the  system  performance 
of  the  adaptive  optics  will  be  limited  by  the  anisoplanatic 
field.  The  calculation  results  in  the  literature  [1,2]  indicate 
that  the  influence  of  the  anisoplanatic  field  on  the  compensation 
effect  is  quite  serious.  To  overcome  the  anisoplanatism  error, 
some  researcher  [3]  proposed  the  concept  of  a  multilayer  adaptive 
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optics  system  in  calibrating  atmospheric  turbulence  in  layers  at 
different  heights.  However,  up  to  now  there  have  been  no  further 
reports  on  theoretical  work  and  experimental  status.  This  paper 
will  further  discuss  the  compensation  for  the  effect  of 
atmospheric  turbulence  by  using  an  adaptive  optics  system  for 
calibration  under  anisoplanatic  conditions;  a  new  method  is  also 
proposed  to  overcome  anisoplanatic  errors. 

II.  Anisoplanatic  Influence  in  Conventional  Adaptive  Optics 

When  the  beacon  deviates  from  the  target,  the  beacon  light 
beam  and  the  emitted  light  beam  will  traverse  two  different 
turbulent  paths.  Strictly  speaking,  even  if  these  two  beams  are 
within  the  same  isoplanatic  field,  the  phase  distortion  of  the 
beacon  channel  and  the  phase  distortion  of  the  emission  channel 
cannot  be  completely  identical.  In  traditional  adaptive  optics, 
conjugated  phase  distortion  in  the  beacon  channel  thus  detected 
is  applied  to  the  initial  emission  field  in  order  to  calibrate 
the  phase  distortion  in  the  emission  channel.  This  compensation 
will  certainly  generate  an  error  in  the  residual  phase. 

Generally  speaking,  the  error  Increases  with  increasing  distance 
between  beacon  and  target*  Assume  that  the  atmospheric 
equivalent  distance  between  beacon  and  target  is  P,  and  the 
target  lies  on  a  symmetrical  axis  of  the  emission  aperture,  the 
effective  thickness  of  the  atmospheric  layer  is  L  (see  Fig.  1). 
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Fig.  1.  Schematic  diagram  of  the  system 
KEY:  1  -  Beacon  2  -  Target  3  -  Turbulent 
atmosphere  4  -  Ground  emission  plane 

Based  on  the  expanded  theory  of  Huygens-Freanel  [4],  the 
average  light  intensity  after  calibration  at  the  target  is 
</<0)>  -2'*L-*jJ/*pl4r,F1«,(F,)a;(F,)exp  </*  > 

•<exp{x(pl.0)  +  Z(/»1.0)  +  /[S(#l.0)-S(pJX))J-a«^F,)-S(^Fi)]}>  <*> 
Here,  u0(p)  is  the  initial-field  distribution;  A  is  the 
wavelength  of  the  emitted  light  beam;  an<*  S(Pi»0)  are, 

respectively,  logarithmic  fluctuation  of  oscillation  amplitude 
and  phase  fluctuation  of  the  light  field  at  the  point  (0,  L) 
after  transmission  through  a  turbulent  atmosphere  for  a 
spherical-surface  wave  at  point  S(P,fi)  is  the  phase 

fluctuation  of  the  spherical-surface  wave  at  point  (1*,P)  arriving 
at  point  (0,fi).  Based  on  the  atmospheric  interchange  principle, 
phase  fluctuation  S(f^,P)  at  point  (L,P)  for  spherical-surface 
waves  of  S(P,fi>  and  the  spherical-surface  wave  at  point  (0,  P*) 
are  the  same.  Both  theory  and  experiments  indicate  that  the 
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fluctuations  In  the  oscillation  amplitude  and  phase  fluctuations 
satisfy  the  Gaussian  distribution;  these  two  fluctuations  are 
independent.  Therefore,  Eq.  (1)  can  be  rewritten  as 

•  exp(  ~  0|(Pi”Pi)  ~  D*t(Pi~Pi  •  ?  )I  (2) 

Here  Dt(pt-p»)  is  the  logarithmic  structure  function  of 
oscillation  amplitude  fluctuations;  f>A.(Pi-Pi«p)  1*  the  structure 
function  of  the  residual  phase  caused  by  anisoplanatism,  and  can 
be  expressed  by  the  following  equation: 

DA,(p,-pi.F)^<[S{pt,0)^S(ppF)-S(pv  0)  +  S(ft.f))*>  (3) 

Function  S(p,P)  can  be  expanded  into  a  Taylor  series.  The 
authors  proved  that  the  higher-order  terms  of  the  expansion 
equation  can  be  neglected  when  satisfying  condition  P  <  0.3D  (D 
is  the  effective  emission  aperture)  by  only  taking  the  single¬ 
order  approximation 

0)  +  ^*  V,S(p,P)|r>4  (4) 

Here,  VfSip,F)  is  the  phase  gradient  along  the  direction  of  ^ 
Thus  the  structure  function  in  Eq.  (3)  can  be  expressed  as 

For  further  computation,  the  authors  expand  a  set  of  orthogonal 
Zernike  polynomials  [5,6]  for  wavefront  phase  distortion 

sip,r)-i»,F,ut/D)  *  (6) 

In  the  equation,  coefficient  aj  reveals  the  proportion  of  mode 
F,(2f/D)  in  phase  distortion.  By  using  Eq.  (6),  we  can  derive 
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(7) 


that  the  phase  gradient  is 

ft 

Substituting  Eq.  (7)  In  Eq.  (5),  if  only  the  first  six  orders 
(manifesting  the  major  function)  are  only  retained  in  the 
expansion  equation,  we  can  obtain 

d»-  y*  tf1 1/, )•?!■)  (it 


With  Kolmogorov's  turbulence  spectrum,  we  have  obtained  the 
correlation  among  various  order  expansion  coefficients;  among 


thesi, 


<ai>-X32xlO-*(-^p  )*• 


Here,  9Q  is  defined  as  the  isoplanatic  angle  of  phase  fluctuation 

9>-[tMk'ld>g**ClU))-”  (l0) 

Substituting  Eq.  (9)  in  Eq.  (8),  we  obtain 


Du-™  Lk  - |£-ft*(I  +  2C') 

rTTi  *V>D'A 


The  above  equation  can  be  rewritten  as 

)“>^-«'(l  +  JC)  (11) 

In  the  equation,  C-eo*<;Mi“£>  J  r0  i#  the  coherent  length  of 

J  •  * 

the  spher lea l«sur face  wave;  9  ■  P/L  is  the  angle  formed  by  the 
beacon  and  target. 

'•■•TST*']*0  T  (,1> 

* 


-*■  •!*■*•  7 — "■ - r  (u) 

Substituting  Eq.  (12)  in  Eq.  (2),  we  obtain 

(IS) 

•  «  «  • 


Assume  that  the  initial  emission  field  is  of  a  gaussian 
distribution 

•».(/») 


(16) 


Here,  F  is  the  radius  of  curvature  of  the  emission  light  beam. 
In  the  following  discussion,  we  will  consider  the  situation 
that  the  light  beam  is  focussed  at  the  target;  that  is,  F  -  Lq. 
Eq.  (15)  can  be  simplified  as 


</(0)>  -  -|jr£r  «P(  -W  >{[*'  expt "  “443( 


W(,  +  2C,)1 


(17) 


In  the  equation,  CTj  is  the  variance  of  the  logarithmic 
fluctuation  of  the  oscillation  amplitude.  After  simple 
integration,  we  can  obtain 

</(0 ))  —  )*(  y\%L%D  ^ 
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Therefore}  we  can  obtain  the  result  that  the  stereo  ratio  Rs 
(characterizing  light  beam  quality)  is 

j»,.{(i+2.2I(  jjh  n[,+6.65(  j. 

<■») 

In  reference  [7],  we  discussed  the  influence  on  the  compensatory 
effect  due  to  fluctuations  of  oscillation  amplitude.  Here  we 
discuss  only  the  effect  of  anisoplanatism.  From  Eq.  (19)  the 
effect  due  to  anisoplanatism  intensifies  with  increase  in  or 
decrease  in  rQ  ;  however,  the  increase  in  aperture  D  can  reduce 
the  effect  of  anisoplanatism.  From  Eq.  (19),  we  obtain  the 
stereo  ratios  of  two  extreme  cases: 


*«-l 


*i-<U2(  A  )*(_L«2.).* 


whm 


Hhen  1 


(20) 

(21) 


From  this  we  can  deduce  that  when  the  included  angle  formed  by 
the  beacon  and  target  is  relatively  large,  and  if  the  aperture  D 
or  the  coherent  length  tq  is  relatively  small,  the  compensation 
for  atmospheric  turbulence  for  the  case  when  conventional 
adaptive  optics  will  be  very  poor;  after  compensation,  light  beam 
quality  may  possibly  be  worse  than  in  the  absence  of 
compensation. 
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III.  New  Method  of  Reducing  Anisoplanatic  Error 


When  there  is  an  included  angle  between  the  beacon  light 
beam  and  the  target  light  beam,  wavefront  distortion  of  the 
beacon  light  beam  cannot  fully  represent  wavefront  distortion 
along  the  emission  path.  Therefore,  the  error  of  the  residual 
phase  will  be  produced  if  one  simply  substitutes  S(P,^))  for 
S(0t( ).  To  reduce  this  anisoplanatic  error,  we  can  indicate  the 
phase  distortion  on  the  target  path  by  using  the  wavefront  phase 

(and  its  gradient)  of  the  beacon  light  beam;  that  is,  by  using 

S-S(P,p)~P' V'S(P.p)  (22) 

as  the  phase  predistortion  applied  to  the  initial  emission  field. 

We  call  this  method  the  phase  gradient  method.  After  the  phase 

calibration  shown  in  Eq.  (22),  we  can  find  that  the  light 

intensity  at  the  axis  is 

•<exp{*(p„  0)  +  x(/>„  0)  +  i[S(p„0)-S(P,p1)  +  ?.  V,S(P,p,)J 
-HS(p1,0)-S(P.pJ)  +  F'  V,J(/\p,)]}>  (23) 

The  above  equation  can  be  simply  rewritten  as 

</  ,(^)'>mX‘1L~l^<PpldipJu%(p,)u*(pi)txp(ik-^^-  ) 

•  exp(  -  ±  Dt(px-pa)  -  \  DU  /»)]  (24) 

In  this  equation,  D£g  is  the  residual-phase  error  under  the  phase 
gradient  method;  the  error  can  be  indicated  by  the  second-order 
derivative  of  the  phase 
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2  ^P‘  dP?  +  2P>P'  dP.dP,  +P'  dP}  ) 


•( S(P ,  Pt)-S(P.  pt)|U)’> 


The  first  ten  terms  are  retained  in  the  phase  expansion  equation, 
then  the  above  equation  becomes 


— ,-14l,<^>  •  !»•  (1  +  4C*) 


At  Kolmogorov's  turbulence  spectrum  <aj> «6.I9x  10*'(  )y*  i 

Lv| 

therefore 

D.>S7.0S(-^  )<•  -Jsr  -<T >“•  l?.-»Ml+4C‘)  (27; 

Substituting  Eq.  (27)  into  Eq.  (24)  and  assuming  that  the 
emission  field  is  the  focused  gaussian  light  beam,  we  can  obtain 


aM)>'W“p< 


-•:)lj<pptxn-2£  -2>.s.(-^  >< 


J5S-  ;<-rfe-)’V- <1+401 


DyVi  '  3.I8L,  '  r  ^  :* 

After  computation  with  integrating  on  the  above  equation,  ve 
obtain 


« '<»)>•  W  «<■<  (I >«( )' ( Tjf^  >M1 


Thus,  we  can  obtain  the  result  that  the  stereo  ratio  in  the  phase 
gradient  method  is 
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,  »*,  ..  3-189  v#  /  Lrt  »Mi«  |  .  ft  <(  *1W  (  *dk  )Mi  )*w 

A,'-{n.+  u-3(-gr-  >’<TTI£d  >  H,  +  7151  0,  ri  3I8V>  '  n 

•exp(-sj)'  (30) 

To  verify  how  effective  the  phase  gradient  method  is,  ve 
should  find  by  how  the  angle  is  expanded  for  effective 
calibration  after  using  the  phase  gradient  method.  With  the 


conventional  method,  the  range  of  angles  for  effective 
calibration  is  the  isoplanatic  angle  9q,  When  the  beam  is  at  the 
fringe  of  the  isoplanatic  angle,  the  stereo  ratio  is 

t  %  •  •  "  ^ 

AM«{(1  +  221  x 3- ),y*l  •(  1  +  6-65 *3.t8vl‘(~^'  YA]  Y'n  (31) 
We  define  the  range  of  angles  for  effective  calibration  in  the 


phase  gradient  method  is  realized  when  the  stereo  ratio  Rg  is 
equal  to  the  9  corresponding  to  Rgo*  Therefore 


A*- U I +M-3  •(-!!$£)•< 


1-lSLfD 

3-189  v.. 


Rs0  is  given  by  in  Eq.  (31).  Therefore,  ve  can  obtain  the  number 
of  times  required  for  expanding  the  angle  range  with  effective 
calibration  after  using  the  phase  gradient  method: 


t 
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In  the  following,  we  consider  the  results  of  homogeneous  and 
inhomogeneous  turbulence  conditions.  For  a  horizontal  path, 
turbulence  is  homogeneous.  In  this  case,  the  atmospheric 
equivalent  distance  P  between  beacon  and  target  is  the  actual 
distance  between  beacon  and  target.  L  is  the  distance  between 
target  and  emission  plane.  According  to  the  definition  of  Lo  , 
we  can  obtain  the  result  that  Lq  ■  3.18L.  Therefore,  on  the 
horizontal  path,  G  is 


( 


f- 


99( 


-g-)4+ 187-03< 


-g-  ),y,l68.74<  -g-  ),/S+  8.86]  -8-48( 


A 

D 


93-52  ( 


Jlvvi 
D  * 


y 

-)  (34) 


On  the  slant  path,  the  turbulence  intensity  along  the  path  is  a 

function  of  height.  Generally  speaking,  turbulence  intensity 

weakens  with  increase  in  height.  As  indicated  in  the 

2 

experimental  results,  Cn  at  a  height  of  20km  is  lower  by  three  or 

four  orders  of  magnitude  than  that  on  the  ground.  Therefore,  the 

major  function  of  anisoplanatism  is  the  atmosphere  below  the 

height  of  201cm.  Above  20km,  6^  basically  does  not  vary. 

Therefore,  when  the  target  is  outside  the  atmospheric  layer  at  a 

height  H  from  ground  surface  (H  >  20km),  we  can  consider  that  the 

effective  thickness  of  the  atmospheric  layer  is  20km.  Therefore, 

if  the  zenith  angle  of  the  target  is  iff,  then  L“20,000  sec  Iff  (m); 

the  atmospheric  effective  distance  P"20,000  •  8  *sec  y’(m)  between 

beacon  and  target;  however,  the  actual  distance  (between  them)  is 

H.0*8ecy.  To  find  the  value  of  Lo»  we  should  know  the 

distribution  of  the  turbulence  intensity  C?  with  height.  Based 

n 
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on  the  accumulated  experimental  data,  reference  [8]  givea  the 
distribution  mode  of  turbulence  height  adaptable  to  different 
situations 

Cl(H)  -A*  H'™*  exp(  -  H/HJ  (  35  ) 

In  the  equation,  Hq  is  the  characteristic  height.  In  the 

-14 

situation  of  weak  turbulence,  A-4«  10  and  HQa2500m.  In  the 

-12 

situation  of  intensive  turbulence,  A"10  and  H0"4OOOm. 
Substituting  Eq.  (35)  into  Eq.  (14)  and  after  computation,  we 
obtain 

£,«3.18sec*.  /Vly(~,  ))-”  (36) 

x  n, 

Here,  y(b,x)  is  an  incomplete  Gamma  function.  When  H  is 
sufficiently  large,  the  above  equation  becomes  a«secy.  In 
the  conditions  of  weak  and  strong  turbulence,  a  is  4428m  end 
7095m,  respectively.  Upon  comparison  with  reference  [5],  these 
values  are  basically  consistent  with  the  values  for  day  and 
night.  Therefore,  we  can  obtain  that  the  amount  of  amplification 
for  the  range  of  angles  in  which  effective  calibration  can  be 
made  under  the  phase  gradient  method  is 


G- 


From  Eqs.  (34)  and  (37),  we  can  see  that  whether  horizontal  or 
slant  path,  the  amount  of  amplification  G  for  increasing  the 
range  of  angles  with  effective  calibration  is  related  to  D/rg  • 
In  Fig.  2,  the  authors  plot  the  calculation  results  of  G  with 
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variation  in  D/rQ.  Obviously,  the  greater  the  D/r q  ,  the  greater 
G  is.  For  identical  D/r q  values,  the  amount  of  amplification  in 
a  horizontal  path  is  larger  than  for  a  slant  path.  For  example, 
when  D-4m  and  r^O.lm,  G  on  the  horizontal  path  is  approximately 
6.6;  however,  on  the  slant  path  corresponding  to  conditions  of 
weak  and  strong  turbulence,  G  is  1.7  and  2.2,  respectively.  We 
should  point  out,  in  addition,  that  when  D/ro  is  sufficiently 
reduced  to  a  certain  value  (under  conditions  of  weak  and  strong 
turbulence,  D/ro  *s  14  and  8.5,  respectively),  the  amount  of 
amplification  G  will  be  smaller  than  1.  In  other  words,  the 
range  of  angles  with  effective  calibration  shrinks  under  the 
phase  gradient  method;  this  is  because  we  only  revised  the  first- 
order  terms  of  phase  gradient  but  neglected  the  higher-order 
terms.  Therefore,  when  D/rQ  is  relatively  small,  the  second- 
order  terms  and  even  higher-order  terms  of  phase  gradient  should 
be  further  considered.  Thus,  the  amount  of  amplification  G  will 
increase;  we  are  carrying  this  research. 
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Fig.  2.  Variation  in  amount  of  amplification 
C  with  D/ro 

horizontal  path _ slant  path 

1  -  Weak  turbulence  2  -  Strong  Turbulence 

IV.  Preliminary  Conclusions 

Anisoplanatism  of  a  beacon  is  an  important  factor 
restricting  adaptive  optics  as  a  technique.  To  effectively 
calibrate  the  effects  due  to  atmospheric  turbulence,  the  beacon 
should  lie  within  the  range  of  the  isoplanatic  angle;  otherwise, 
the  compensatory  effect  will  be  degraded.  Generally,  the 
isoplanatic  beam  is  very  small  (only  about  several  seconds  in 
angular  measure);  this  imposes  rigorous  requirements  on  the 
beacon.  The  authors  conducted  some  improvements  on  aspects  of 
the  principle  to  significantly  increase  the  range  of  angles  for 
effective  calibration,  thus  the  effect  of  anisoplanatic  error  is 


greatly  reduced.  As  shown  by  the  analytical  results:  along  the 
horizontal  path,  the  effective  calibration  range  of  several 
seconds  of  original  angular  measure  is  expanded  to  tens  of 
seconds  by  using  the  phase  gradient  method  with  4m  aperture. 

Thus,  the  anisoplanatic  problem  of  the  beacon  is  solved.  Along 
the  slant  path,  the  range  of  angles  with  effective  calibration 
can  also  be  increased  by  about  twofold.  By  using  an  even  larger 
aperture,  the  problem  of  beacon  anisoplanatism  can  be  solved.  In 
particular,  in  the  infrared  waveband,  the  problem  can  be 
basically  solved  by  a  factor  of  2  taken  as  the  amount  of 
amplification.  In  practical  execution,  only  with  a  beacon  light 
source  (such  as  a  moving  target  or  multistar  system)  can  we 
provide  for  phase  gradient  of  real  time  measurements,  so  it  ik 
not  difficult  to  carry  out  the  calibration  v  der  the  phase 
gradient  method. 

In  writing  this  paper,  the  authors  had  profitable  discussion 
with  Jiang  Uenhan,  and  express  gratitude  to  him. 

The  first  draft  of  the  paper  was  received  on  27  July  1990; 
the  final  revised  draft  was  received  for  publication  on  14 
September  1990. 
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